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ABSTRACT 

A  new  Generalized  Cross-Correlator  (G^CjCj)  for  the  passive 


time  delay  estimation  problem  is  presented.  The  interpreta¬ 
tion  of  this  G,C,C.  is  that  of  estimating  the  cross-correlation 
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function  by  cross-correlating  the  least  mean  square  esti¬ 
mates  of  the  signal  component  in  each  of  the  observed 
waveforms.  The  implementation  is  simply  a^C,  with  the 
weighting  Alter  equal  to  the  magnitude  coherency  squared. 
Numerical  evaluation  of  the  performance  of  this  processor 
and  a  robust  version  indicates  that  they  compare  favorably 
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INTRODUCTION 


The  estimation  of  propagation  delay  in  a  common  signal  arriving  at 
two  spatially  separated  sensors  is  a  problem  which  has  received  much 
attention  in  the  literature.  Cross-correlation  methods  are  particularly 
popular  because  of  the  richness  and  variety  of  processors  within  this 
class  and  the  general  ease  of  implementation.  Of  these,  the 
Hannan  Thompson  {H.T.)  [1,5],  the  Eckaert  [l]  and  the 

Hassab -Boucher  ( H.B .)  [2]  are  examples  of  "optimum"  processors 
which  maximize  some  performance  criteria.  On  the  other  hand,  the 
SCOT  [7]  and  the  simple  cross-correlator  (C.C. )  are  examples  of  "Ad  Hoc" 
or  intuitive  correlation  type  processors.  In  general  the  optimum  proces¬ 
sors  are  very  sensitive  to  deviations  from  the  assumed  signal  and  noise 
characteristics.  By  way  of  contrast,  the  C.C.  and  SCOT  appear  to  be 
more  robust  to  these  deviations  from  the  nominal.  However,  these  later 
processors  can  have  very  poor  performance  at  the  nominal  point. 

In  the  following  development  an  intuitive  approach  to  the  problem 
yields  another  type  of  (generalized)  cross-correlator.  This  is  the 
Wiener  Processor,  or  W.P. ,  which  has  a  simple  form.  Yet  preliminary 
results  indicate  it  outperforms  most  of  the  other  above  named  proces¬ 
sors  when  compared  under  various  performance  criteria  for  the  few 
important  cases  considered  in  this  paper.  A  "robust"  version  of  the  W.P. 
also  indicates  potentially  good  performance  relative  to  the  others  under 
uncertain  spectra. 


L  PROBLEM  STATEMENT  AND  BACKGROUND 


We  first  consider  a  system  model  generating  the  observations  in  Fig. 
I.  We  observe  Gaussian,  ergodic,  wide-sense  stationary  processes  xx (t) 
and  xz(t)  over  a  time  interval  [0,T]  which  contain  uncorrelated  noises 
nj(0  and  n^t)  and  signals  s(t )  and  s0(f)  respectively.  We  assume  c(t) 
is  a  linear  time  invariant  channel  having  a  transfer  function  C(cj)  with 
unknown  linear  phase  so  that  s9(t)  is  a  delayed  but  possibly  distorted 
version  of  s(f ).  Furthermore,  we  assume  that  the  noises  are  uncorre¬ 
lated  with  the  signal  and  that  T  is  much  greater  than  the  correlation  time 
of  Xj(£)  and  x2(f  )•  The  object  is  then  to  estimate  the  time  delay,  D.  asoci- 
ated  with  the  channel. 

We  define  the  sample  cross-correlation: 

ki2(T)=  ^f£iz{u)ejvTdu  .  (1-0 

Here  &t2(o)  is  some  unbiased,  consistent  estimate  of  the  cross-spectrum 
Ciziu).  For  definiteness  we  assume  that  £«(«)  is  obtained  by  the 
Bartlett  Procedure  [19]  consisting  of  segment  averaging  periodogram 

type  estimates  ^-ArJk(«)A'ak(cj)  where  m  is  the  total  number  of  subdivi¬ 
sions  of  [O.T]  and  *i*(u)  and  X&iu)  are  the  Fourier  Transforms  of  *j(f) 
and  x2(f )  on  the  fc-th  segment. 

For  very  large  observation  time  and  suitable  m,  ftjg(r)  gives  a  good 
approximation  to  the  true  cross-correlation  function  which  has  a  global 
peak  at  D.  In  fact  if  c(f)  is  pure  delay  and  #(f)  is  white,  the  cross- 
correlation  function  is  a  delta  function  at  the  true  delay.  For  finite 


observation  time  we  can  decompose  ft12(r)  into  the  sum  of  four  terms: 


fttffr)  =  c(t)*&-(t)+  c(t)*^(t)  +  ft^r)  +  (1.2) 

where  denotes  convolution.  Here  (t)  is  an  estimate  of  the  signal 
autocorrelation  function  R„(t)  and  ^.(r),  ft^r)  and  ^^(t)  are  esti¬ 
mates  of  the  cross-correlation  between  the  signal  and  noise  terms  in  the 

observations.  In  the  limit  the  sample  cross-correlation  converges  to 

♦ 

c  (r)  *Rm(r)  which  displays  an  absolute  maximum  at  D.  Thus,  it  is  the  last 
three  terms  in  Eqn.  (1.2)  which  constitute  zero  mean  disturbances 
affecting  peak  resolution  of  the  first  term.  This  suggests  prefiltering  the 
sample  cross-correlation  with  a  filter  W(cj)  to  obtain  better  resolution  of 
the  peak  at  D,  where  JP(cj)  has  zero  phase.  This  scheme  is  referred  to  as 
the  generalized  cross-correlation  method  or  the  Generalized  Cross- 
Correlator  (G.C.C.)  and  is  illustrated  in  Fig.  II  (The  content  of  the  dashed 
box  is  only  symbolic  for  the  operation  in  Eqn  (1.1)).  We  denote  the  G.C.C. 
output  waveform  Rfc(r).  Therefore  we  have: 

Rfe(r)  =  J^fZK(»m»)'iurdU  .  (1.3) 

When  W(u)  is  unity  the  resulting  G.C.C.  is  called  the  simple  cross- 
correlator  (C.C.).  Considering  the  first  term  in  Eqn.  (1.2)  as  a  ''signal”  in 
additive  noise,  classical  optimal  filtering  theory  can  be  applied  to  derive 
filters  W(cj)  which  maximize  signal-to-noise  ratio. 

Letting  the  last  three  terms  of  Eqn  (1.2)  be  characterized  as  "noise” 
we  can  define  a  signal-to-noise  ratio  as  the  energy  at  the  output  due  to 


-  4- 

the  global  peak  of  c  (r)  'ft*  (t)  divided  by  the  output  power  due  to  the 
three  noise  terms.  We  will  denote  this  by  SNR1.  For  the  specific  cross- 
spectral  estimation  procedure  outlined  above,  the  noise  power  has  been 
derived  and  is  given  by: 

(1.4) 

k  is  a  positive  constant  inversely  proportional  to  T,  the  observation  time. 
Gn(o)  and  Gzku)  are  the  power  spectral  densities  of  Ihe  observations 
Xi(t )  and  x2(t)  respectively.  7?2(“)  is  the  magnitude  coherency  squared, 
where,  defining  the  true  cross-spectrum  G^tS), 


7?2(w)  = 


\Gi2(u)\* 

Gu  {o)Gzt(u) 


Then  from  the  defining  relation: 


(15) 


SNR i  = 

we  obtain: 

(16) 

m 

[^-/  \Gizmw(U)dU)z 

SNR ■,  =  -2^= - 

gjj-/  C11(u)Ca2(cj)(l-7f2(w))|W(«)|2d« 


Hie  maximum  is  obtained  through  the  Schwarz  inequality  and  yields  the 
H.  T.  processor.  The  same  result  is  derived  in  [5]  as  the  result  of  minim¬ 
izing  the  local  variance  of  the  delay  estimate  over  the  entire  G.C.C.  class, 
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and  in  [1]  as  the  result  of  maximum  likelihood  estimation.  The  filter  is: 


*U.T.(u) 


1  7\ziu) 
\Giz(u)\  1  —y  12(d) 


(IT) 


Neglecting  the  effect  of  the  signal  and  noise  cross  terms  , 
c(T)#ft**(T)  and  R^u^t)  in  Eqn.  (1.2),  gives  another  characterization  of 
the  noise  in  the  cross-correlation  domain.  With  this  definition  of  noise 
another  signal-to-noise  ratio  is  defined  in  [8],  SNR 2,  which  is  shown  to  be 
maximized  by  the  Eckart  processor  fVg^(cj): 


m 

t£rf  \Gn{u)\W{u)du? 

SNR  z  =  - =-^= -  (1.8) 


Q»l(w)G>»«(w) 


(19) 


Where  £*,(&>)  and  G^cj)  are  the  auto-spectra  of  the  noises  nj(f)  and 
n^t)  respectively.  Note  in  terms  of  the  spectra  of  the  observables  *j(f) 
and  x2(0  the  filter  takes  the  form 


m  (f.\  _  _ _ 

a  (Cn(w)  -  |Cia(w)|)(Cj52(w)  -  |Cj2(Q)i) 

which  attests  to  the  complexity  of  implementing  the  Eckart  in  general. 

Hassab  and  Boucher  [2]  take  the  approach  of  maximizing  a  signal-to- 
noise  ratio,  SNR*  defined  as  the  ratio  of  the  expected  peak  energy  at  the 
true  delay  to  the  total  statistical  variation  of  the  output  of  the  C.C.C. 
This,  in  some  sense,  lumps  the  "signal",  c  (t)  'ft,,  (t),  variation  into  the 


-  0  - 

noise  terms  and  yields  the  H.B.  filter  WH  p  (&>): 


l-£rf  |Ci2(u)|H,(w)dw]8 

SNR s  =  -4" -  (110) 

|r/c11(«)G28(«)  \W(u)\zda 


*U.B.  (a)  = 


|Glg(M)| 

Gn{cS)Ggz(u) 


(111) 


The  H.B.  is  similar  to  the  SCOT  introduced  by  Carter  si  a l  [7]  in  that  it 
suppresses  the  cross-spectral  estimate  in  u- regions  of  both  high  and  low 
signal-to-noise  ratio  in  an  attempt  to  reject  strong  tonals  in  the  observa¬ 
tions. 


r 
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H.  THE  WIENER  PROCESSOR 

Here  a  slightly  different  approach  is  taken  to  derive  an  optimal  filter. 
We  deal  directly  with  the  quantities  in  the  observation  time  domain  (i.e. 
Fig.  I).  The  procedure  is  motivated  by  the  following  argument.  If  we 
knew  the  signal  s(f)  and  the  filtered  version  s,  (0  exactly  then  from  the 
linearity  of  the  phase  of  the  channel  the  time  delay  could  be  estimated 
exactly  by  detecting  the  peak  of  the  sample  cross-correlation  of  s  (f )  and 
s0  (0-  Therefore  we  simply  try  to  estimate  the  signal  s[t)  as  best  we  can 
from  the  observations  Xj(f )  and  the  channel  output  signal  ss(£)  from 
x2(t)  by  minimizing  the  mean  square  errors: 

£'{(s(£)  -  s(0)2$  =  min  (2.1) 

El(so  (0  =  min  (2.2) 

where 

r 

s(£)  =  f  xi{a)hl{t'-  a)da 
-T 
T 

»«(£)  =  f  x2{o)h2(t-  cr)d<7  . 

-T 

The  above  procedure  is  illustrated  in  Fig.  III.  Given  the  channel  charac¬ 
teristic  C(cj)  the  solutions  to  Eqn.  (2.1)  and  (2.2)  are  the  Wiener  filters 
Hi(u)  and  H2(u): 


= 


<*•(«)  + 


//2(«)  = 


C»fcK(»)i8 

C„(U)|C(U)|2+  ^.(cj)  * 


Noting  that  G]2(u)  =  C(u)  £*,(&>)  we  can  egress  the  above  filters  in  terms 
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of  the  quantities  derived  from  the  observables: 


H1(u) 


1 

C(u) 


H2(u)  =  C*(u ) 


Gizfa) 

Gn(«) 

Giz(u) 

^22(0) 


(2.3) 

(24) 


where  C*(u)  is  the  complex  conjugate  of  C(u). 

With  these  filters  the  sample  cross-correlation  of  the  least  mean 
square  error  estimates  of  s  (£ )  and  se(f)  yields  the  estimate  of  the  cross 
correlation  function: 


(2.5) 

R£p(t)  =  jr|  ^(cj)S0k(U)e^dU 

where 

5*(<y)  =  /f1(«)A'u(  u)  (2.6) 

S,4(«)  =  (2.7) 

Xuciu)  and  A'at(cj)  are,  as  before,  the  finite  time  Fourier  transforms  of 
*2 (t)  and  i2(0  on  the  fc-th  segment  of  [0,T].  Regrouping  terms  in  (2.5) 
we  obtain 


=  £-/  6-2<“)  eS"du  (2-8) 

Comparing  Eqn.  (2.8)  with  Eqn.  (1.3)  we  have  the  result  that  the  W.P.  is 
equivalent  to  using  a  Generalized  Cross-Correlator  with  the  filter  W  (cj) 
equal  to  the  magnitude  coherency  squared. 

It  should  be  emphasized  that  even  though  the  Wiener  filters  Hi  and 
Hz  involve  the  knowledge  of  the  channel  C(u)  itself  the  G.C.C.  equivalent 
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processor  does  not  impose  this  requirement.  In  fact,  as  far  as  the  cross¬ 
correlation  estimate  of  time  delay  is  concerned,  the  actual  channel  is 
immaterial  to  the  peak  detection  procedure  in  the  cross-correlation 
domain.  Hence  the  Wiener  filter  implementation  (Fig.  Ill)  with  C(u )  arbi¬ 
trarily  set  to  unity  in  Eqns.  (2.3)  and  (2.4)  is  equivalent  to  any  other 
choice  of  C(cj)  for  the  time  delay  estimation  problem. 

The  definition  of  "additive  noise"  leading  to  the  signal-to-noise  ratio 
SNR j,  Eqn.  (1.4),  yields  the  interpretation  that  yfz (w)  is  a  measure  of  the 
proportion  of  the  power  spectral  density  of  the  sample  cross-correlation, 
Rj2(t),  which  is  due  to  the  sample  autocorrelation  ^(r)  (the  first  term 
in  Eqn.  (1.2)).  Thus  the  W.P.  de-emphasizes  those  cj-regions  where  the 
sample  cross-spectrum  is  most  likely  to  be  a  highly  inaccurate  estimate 
of  the  true  cross-spectrum.  This  is  not  surprising  given  the  raison  d'etre 
of  the  W.P.  which  is  to  accurately  estimate  the  smoothed  sample  auto¬ 
correlation,  R*,(t)*c  (t). 

The  W.P.  does  not  of  course  maximize  the  signal-to-noise  ratio  in 
general.  If  we  examine  the  optimal  processor  for  SNR  j,  the  H.T.  (Eqn. 
(1.7)),  we  see  that  it  has  the  additional  ability  to  overemphasize  as  well  as 
to  de-emphasize  the  cross-spectral  estimate  according  to  the  function 
7?2(o)/(1  -7?2(cj))-  (Actually  in  [5]  the  above  function  is  shown  to  be 
inversely  proportional  to  the  variance  of  the  phase  estimate 
&i2(«)/|&i2(«)!  with  respect  to  the  true  phase  of  the  cross-spectrum). 
However,  in  situations  where  the  coherence  is  low,  and  the  signal  spec¬ 
trum  nearly  flat,  the  H.T.  and  the  W.P.  are  virtually  identical  and  exhibit 
identical  performance  (Eqn.  (1.7)  becomes  proportional  to  yf2). 


- 


It  is  also  observed  that  the  W.P.  is  equivalent  to  the  H.B.  for  nearly 
flat  signal  spectra  and  also  to  the  Eckart  if  we  add  a  low  signal-to-noise 


4 

* 
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m.  PRACTICAL  CONSIDERATIONS 

In  practice  the  various  spectral  quantities  that  appear  in  the  expres¬ 
sions  for  the  optimal  Alters  are  unknown  (either  partially  or  totally)  so 
that  the  Alters  are  not  directly  implementable.  Two  approaches  to  this 
problem  are  possible.  We  either  estimate  the  spectra  and  substitute  the 
estimates  into  the  aforementioned  Alters  (totally  unknown  spectra)  or  we 
search  for  a  robust  solution  over  a  range  of  spectra  perturbed  from  some 
nominal  point  (partially  unknown  spectra).  In  this  section  the  estimation 
approach  is  brieAy  discussed  and  in  the  next  section  the  robust  approach 
is  considered. 

Several  simulation  studies  have  recently  appeared  in  the  literature 
which  exercise  the  various  Alters  for  known  signal  and  noise  spectra.  To 
the  extent  that  generalizations  are  possible  from  a  limited  set  of  simula¬ 
tions.  they  have  shown  that  the  H.  T.  processor  seems  to  outperform  the 
others  such  as:  the  SCOT  and  simple  cross-correlator  [10];  the  Eckart] 
and  the  H.B.  for  moderately  broadband  signals  [11,12], 

When  the  spectral  quantities  are  totally  unknown  the  H.T.  is  no 
longer  implementable.  Two  main  approaches  to  optimal  time  delay  esti¬ 
mation  have  been  investigated.  The  exact  maximum  likelihood  procedure 
for  unknown  spectra  was  derived  in  [9]  by  Hannan  and  Thompson  but  it 
involves  the  simultaneous  solution  of  several  coupled  equations  at  a  large 
number  of  grid  points  in  frequency.  The  other  approach,  the  Haman 
Hannan  procedure  (H.H.)  [4],  directly  substitutes  spectral  estimates 
into  the  H.T.  to  obtain  the  G.C.C.  structure: 
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= tZ  ,3So  *f0<“>  *  "T,d“  (3, 0 

where  7f2(cj)  is  an  estimate  of  magnitude  squared  coherence  and  is 
an  estimate  of  the  phase  ?(cj)  of  the  cross-spectrum  between  the  obser¬ 
vations.  This  procedure  converges  to  the  H.  T.  in  the  limit  of  large  obser¬ 
vation  times  for  suitable  estimates  of  yfz  and  <p. 

The  H.H.  procedure  was  simulated  in  [10]  and  for  the  specific  spec¬ 
tral  estimates  reported,  the  H.H.  (therein  referred  to  as  the  AM.L. )  was 
shown  to  perform  badly  even  with  respect  to  the  (C.C.)  for  lowpass  and 
bandpass  flat  signed  spectra.  Indeed  a  simple  local  analysis  of  the  filter 
estimation  error  for  the  H.H.  filter,  (cj)  =  3?2 (u)/(l- 7?2(«)),  yields 
the  variance: 

w ( WH.H.  (»))  "  '  VOr(7?2(u))  • 

For  the  specific  coherence  estimate  developed  in  [13]  by  Carter  and  used 
in  Carter’s  program  for  time-delay  estimation  [  14]  the  filter  estimate  has 
a  variance: 

/«,  ,  yhM 

where  Kx  is  a  positive  constant.  Therefore  the  H.H.  may  vastly  under¬ 
weight  those  u-regions  where  coherence  is  highest  because  of  the  large 
slope  of  ITjfif.  as  a  function  of  the  coherence  estimate  when  7?z  is  near 
unity.  In  addition,  using  averaged  type  coherence  estimates  like  that  in 
[13]  ?f2(o)  can  exhibit  a  large  negative  bias  when  the  phase  of  GJ2(«)  is 
rapidly  varying  [4].  This  additional  bias  term  in  Eqn.  (3.2)  could  lead  to  a 
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drastic  increase  in  the  error  of  the  filter  estimate  in  w-regions  of  high 
coherence.  The  net  effect  would  be  to  decrease  output  signal-to-noise 
ratio,  dramatically  raising  the  variance  of  the  time  delay  estimate. 

The  W.P.  on  the  other  hand  can  be  implemented  by  directly  filtering 
the  cross-spectral  estimate  with  the  coherence  estimate,  the  resultant 
filter  having  potentially  less  bias  in  both  the  low  and  high  coherence  u~ 
regions  than  the  H.H.  filter.  This  improved  filter  estimate  could 
translate  into  improved  time  delay  estimator  performance  in  some  cases. 
The  practical  performance  of  this  implementation  of  the  W.P.  remains  to 
be  investigated. 
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IV.  A  ROBUST  VERSON  OF  THE  WIENER  PROCESSOR 

When  the  spectra  underlying  the  observations  are  only  partially 
known  a  rich  theory  developed  in  the  last  decade  [15-17]  can  be  used  to 
design  robust  Alters  for  classical  matched  or  Wiener  Altering.  These  are 
maximin  Alters,  i.e.  they  achieve  the  best  possible  performance  over  all 
spectra  in  a  given  class.  Usually,  one  Arst  Ands  the  signal  and  noise  spec¬ 
tral  pair  which  gives  the  least  favorable  performance,  if  it  exists.  Then 
one  optimizes  the  Alter  for  the  least  favorable  point,  hence  the  name 
maximin  Alter. 

For  the  time  delay  problem  the  only  published  result  known  to  us  in 
this  category  of  processors  is  that  of  Kassam  and  Hussaini  [18]  where 
they  used  the  fact  that  the  Eckart  processor  maximizes  a  classically 
deAned  signal-to-noise  ratio  to  relate  the  Altering  problem  to  robust 
hypothesis  testing.  This  is  achievable  only  by  associating  uncertainty 
classes  with  the  spectral  product  €*,(&>)•  rather  than  with  the  indi¬ 
vidual  noise  spectra  themselves.  Finding  robust  solutions  to  the  other 
signal-to-noise  ratios  SNRZ,  SNRs,  remains  an  only  partially  solved  prob¬ 
lem. 

For  the  pure  delay  channel  an  alternate  approach  to  combatting 
against  poor  performance  with  uncertain  spectra  is  suggested  by  the 
recent  work  in  robust  Wiener  Altering  [18,  17]  when  applied  to  the  W.P. 
With  regard  to  the  original  formulation  of  the  H  P.  we  can  replace  the 
least  mean  square  estimates  of  the  channel  input  and  the  channel  output 
by  the  robust  least  mean  square  estimates  of  s(f )  and  st{t)  under 
uncertainty  in  the  signal  and  noises.  SpeciAcally  we  assume  that  the 
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signal  spectrum  Gj*(cj)  belongs  to  the  spectral  class  joj,  and  that  the 
noise  spectra  £*,(0)  and  G^f(cj)  belong  to  the  spectral  classes  (t)i|  and 
{772}  respectively.  Then  we  solve  for  the  least  favorable  pairs  for  Wiener 
Altering  |£*^(u)  ,  (cj){  and  (<^(cj)  ,  over  the  product  classes 

iaxrj jJ  and  loxrjzi  which  yield  the  robust  Wiener  filters  H f  and  H§  (see 
Eqns.  (2.3)  and  (2.4)).  Finally  we  implement  these  filters  in  the  cross- 
correlation  domain  as  a  G.C.C.,  a  scheme  which  we  will  call  the 
Robust  Wiener  Processor  or  the  R.  W.P, 


-  16- 


V.  NUMERICAL  COMPARISONS 

At  the  present  time  no  simulation  results  concerning  the  experimen¬ 
tal  performance  of  the  W.P.  and  R.W.P.  as  opposed  to  the  other  G.C.C.’s 
are  available.  In  their  absence  a  preliminary  investigation  of  the  relative 
merits  of  the  above  processors  was  performed  based  on  various  signal- 
to-noise  ratio  criteria  for  some  specific  observation  spectra  and  for  the 
pure  delay  channel. 

Fig.  IV  -  Vin  show  the  relative  performance  of  the  H.  T.  ,H.B.  .Eckart, 
SCOT,  and  C.C.  under  the  criteria  SNRlt  SNRz,  SNRz  and  local  variance 
of  the  time  delay  estimate  [3],  for  a  third  order  Markov  signal  in  first 
order  Markov  noises  with  the  noise  3dB  bandwidth  a  factor  of  ten  greater 
than  that  of  the  signal.  The  interesting  thing  to  note  is  that  under  SNR1 
and  SNR 2  the  W.P.  exhibits  better  performance  than  all  of  the  other 
sub-optimum  G.C.C.'s  for  that  particular  definition  of  SNR  and  under 
SNRz  is  in  for  a  close  second  next  to  the  li.L.E.  In  fact,  under  the  cri¬ 
terion  SNRx  performance  of  the  W.P.  is  virtually  identical  to  the  the 
optimal  H  T.  processor.  The  local  variance,  although  it  ranks  the  W.P. 
behind  the  H.T.,  H.B.,  and  Eckart  processors,  only  marginally  disfavors 
the  W.P  at  low  signal-to-noise  ratios.  (It  is  to  be  noted  from  Fig.  VIII  that 
the  SCOT  and  the  C.C.  have  local  variance  orders  of  magnitude  worse 
than  the  W.P  ). 

In  Figs.  DC  and  X  the  performance  of  the  R.W.P. ,  W.P.  and  other 
G.C.C.’s  are  compared  using  SNR j  for  the  e-contaminated  uncertainty 
class  on  the  specific  spectra  in  the  example  outlined  in  Kassam  and  Urn’s 
paper  on  Robust  Wiener  filtering  [16].  Specifically,  under  the  nominal 
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assuxnption  at  each  sensor  we  have  a  signal  with  the  flat  bandlimited 
spectrum  G&(o)  in  first  order  Markov  noise  with  the  spectrum  Gg(u), 
where  the  signal  and  noises  are  of  comparable  bandwidths.  The  uncer¬ 
tainty  on  the  signal  and  noise  spectra  are  modeled  as  the  e-mixtures 

(1-Ci){4(m)  +  ej(i(u) 

and 

(1-c2)GS;(cj)  +  c2Gh(w) 

respectively  with  {*,(«)  and  G^(cj)  arbitrary  spectra  having  the  same 
mass  as  the  nominal  and  et  and  e2  lying  in  the  interval  [0,1].  Fig.  DC 
shows  the  relative  performance  for  the  nominal  spectra  and  Fig.  X  the 
performance  for  the  least  favorable  signal  and  noise  spectra  for  Wiener 
filtering  when  =  0.2  and  e2  =  0.1.  Looking  at  the  nominal  case  (Fig.  DQ. 
we  note  that  the  use  of  the  R.  W.P.  entails  a  loss  of  about  3dB  at  low  SNR 
(below  about  Odb)  over  the  optimal  for  the  least  favorable  pair.  However 
when  the  true  signal  and  noise  spectra  are  least  favorable  for  Wiener 
filtering  the  R.W  P,  displays  uniformly  better  relative  performance  gain¬ 
ing  about  3db  over  the  other  processors  at  low  signai-to-noise  ratios. 
Note  that  this  is  not  necessarily  the  least  favorable  pair  for  H  T  filtering 
so  that  no  conclusive  result  is  indicated  here  However,  Fig  X  does  sug¬ 
gest  that  at  least  for  some  spectra  in  the  al  '>ve  uncertainty  class  we  can 
expect  better  performance  with  the  R  W.  P.  than  with  the  optimal  scheme 
for  the  nominal  spectra. 
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VI.  CONCLUSION 

We  have  outlined  the  development  of  several  of  the  most  popular 
G.C.C.  implementations  and  have  introduced  another  G.C.C.  which  is 
believed  to  be  new.  This  G.C.C.  should  be  easier  to  implement  than  all 
but  the  simple  cross-correlator,  because  the  frequency  weighting  can  be 
estimated  directly  for  unknown  spectra.  Classical  robustness  theory  led 
to  a  simple  alteration  of  the  W.P. ,  which  we  called  the  R.W.P..  The 
evaluation  of  theoretical  performance  in  several  specific  signal  and  noise 
environments  was  presented  which  suggests  that  the  W.P.  and  R.W.P. 
may  be  viable  alternatives  to  existing  time  delay  estimation  schemes. 
Further  experimental  and  simulation  based  performance  evaluation  is 
required  before  any  general  conclusions  can  be  drawn. 
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